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The sector decomposition approach to real radiation at NNLO 
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A method based on sector decomposition is briefly described which is suitable for the calculation of the double 
real radiation part of e + e~ — > 3 jets at next-to-next-to-leading order in a fully differential way. 
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1. Introduction 

The successful comparison with theory of a 
wealth of data collected at high energy collider 
experiments in the past decades largely relied on 
the calculation of radiative corrections. Calcu- 
lating at next-to-leading order (NLO) accuracy 
in perturbation theory is in many cases sufficient 
to match the experimental precision. However, 
there are also prominent exceptions where the 
NLO predictions still suffer from large uncertain- 
ties, or where the experimental precision is ex- 
tremely high. The latter is the case for mea- 
surements of jet rates and shape observables in 
e + e _ annihilation, which for example allow for a 
precise determination of the strong coupling con- 
stant a s . As only NLO predictions are available 
for e + e~ — > 3 jets, the corresponding LEP mea- 
surements have not been included in the current 
world average for a s . A future International Lin- 
ear Collider (ILC) will allow for precision mea- 
surements at the per-mille level, which offer the 
possibility of a determination of a s with unprece- 
dented precision. However, this will only be pos- 
sible if theoretical predictions for jet rates and 
related observables at NNLO accuracy are avail- 
able. 

The calculation of e + e~ — > 3jets at order aj? 
requires the calculation of virtual two-loop correc- 
tions combined with a 1 — > 3 parton phase space, 
one-loop corrections combined with a 1 — > 4 par- 
ton phase space where one parton can become 
soft and/or collinear ("theoretically unresolved"), 
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and the tree level matrix element squared for 
1 — ► 5 partons where up to two partons can 
become unresolved. After the virtual two-loop 
corrections have become available yQ , the stum- 
bling block now is given by the double real radi- 
ation part, where the unresolved particles lead 
to a complicated infrared singularity structure 
which manifests itself upon phase space integra- 
tion. These singularities have to be subtracted 
and cancelled with the ones from the virtual con- 
tributions before a Monte Carlo program can be 
constructed. To achieve this task, two different 
approaches have been followed, one relying on the 
manual construction of an analytical subtraction 
scheme |2I3I4I5I6I7I8I9| , the other one relying on 
sector decomposition |1 ()I1 111 211 311411 .511 till VII Sj 
The application of sector decomposition [191201 
12] to NNLO phase space integrals has first been 
proposed in |10| . Subsequently, it has been ap- 
plied to calculate inclusive phase space integrals 
for e+e- -> 2jets at NNLO |llll3ll4| . The com- 
bination of the sector decomposition approach 
with a measurement function first has been pre- 
sented in |12j . and already lead to a number of 
important results |15I16I17| . Its application to 
the double real radiation part of e + e~ — > 3jets 
at order oP s ^B] represents a new degree of com- 
plexity and therefore will be a crucial test of the 
potential of this method. 

The main features of the two approaches are 
the following: Within the methods based on the 
explicit construction of a subtraction scheme, the 
subtraction terms are integrated analytically over 
the unresolved phase space such that the pole co- 
efficients are obtained in analytic form. This re- 
quires appropriate phase space factorisation and 
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subtraction terms which are simple enough to be 
integrated analytically in D = 4 — 2e dimensions. 
The method naturally leads to a close to minimal 
number of subtraction terms, and observing the 
way how the poles of the different contributions 
cancel allows insights into the infrared structure 
of QCD. 

In the sector decomposition approach, the poles 
are isolated by an automated routine and the pole 
coefficients are integrated numerically. The ad- 
vantages of this approach reside in the fact that 
the extraction of the infrared poles is algorith- 
mic, and that the subtraction terms can be ar- 
bitrarily complicated as they are integrated only 
numerically. On the other hand, the algorithm 
which isolates the poles increases the number of 
original functions and in general does not lead to 
the minimal number of subtraction terms, thus 
producing rather large expressions. 

2. The Method 

The universal applicability of sector decompo- 
sition goes back to the fact that it acts in param- 
eter space by a simple mechanism. The param- 
eters can be Feynman parameters in the case of 
multi-loop integrals, or phase space integration 
variables, or a combination of both. In the fol- 
lowing, the working mechanism of sector decom- 
position will be outlined only briefly, details can 
be found in |2HT3] . 

An overlapping singularity in parameter space 
is of the type 



I = 



dx dy x 1 e (x + y) 1 , 
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where a naive subtraction of the singularity for 
x — > fails. To solve this problem, one can split 
the integration region into sectors where the vari- 
ables x and y are ordered by multiplying Q with 
unity in the form [0(x — y) + Q(y — x)]. Then 

(a) (b) 

the integration domain is remapped to the unit 
cube: after the substitutions y = x t in sector (a) 
and x — yt in sector (b), one has 



dx x 



dyy~ 



dtr 



H + t)-\ (2) 



where the singularities are now factorised. For 
more complicated functions, several iterations 
of this procedure may be necessary, but it is 
easily implemented into an automated subrou- 
tine. Once all singularities are factored out, they 
can be subtracted and the result can subsequently 
be expanded in e. Note that the subtractions 
of the pole terms naturally lead to plus distri- 
butions by the identity 
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where 
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dxf(x) [g(x)/x} + = fdx I[x) fi0) g(x) . 
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In this way, a Laurent series in e is obtained, 
where the pole coefficients are sums of finite pa- 
rameter integrals which can be evaluated nume- 
rically. 

For the numerical evaluation of loop integrals 
it has to be assured that no integrable singulari- 
ties are crossed which spoil the numerical conver- 
gence. For integrals depending only on a single 
scale, which can be factored out, this does not 
pose a problem at all. For integrals with more 
than one scale, like for example two-loop box dia- 
grams, the situation is more difficult, but in the 
case of e + e~ annihilation to massless final state 
particles, evaluation over the whole physical re- 
gion is possible, as the kinematics of these pro- 
cesses is such that the Mandelstam variables are 
always non-negative. 

3. Application to e+e -> 3 jets at NNLO 

In order to focus on a concrete example of phe- 
nomenological relevance, we will discuss the ap- 
plication of sector decomposition to the calcula- 
tion of e + e~ — > 3jets at NNLO in the follow- 
ing. Although the virtual contributions are not 
the main subject of this talk, they will be com- 
mented on briefly. 
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3.1. Double virtual and mixed real- virtual 
contributions 

The contributions to the amplitude which in- 
volve virtual integrals are composed of the two- 
loop (and one-loop times one-loop) corrections 
combined with a 1 — > 3 particle phase space, and 
the mixed real-virtual contributions where one- 
loop corrections are combined with a 1 — > 4 par- 
ticle phase space with up to one unresolved parti- 
cle. In both cases, sector decomposition for loop 
integrals |2J can serve to extract the poles in 1/e. 

In the case of the two-loop integrals, the re- 
sult will only depend on the invariants y± = 
Si2/q 2 ,U2 = si 3 /q 2 and y 3 = s 23 /q 2 , where q 2 
is the invariant mass of the e + e~ system and 
Ei=i V% = !■ However, as the full two- loop ma- 
trix element is known analytically 1 , this part 
could also be taken from the literature, thus gain- 
ing a considerable amount of CPU time. The 
subsequent phase space integration over the yi is 
trivial, and the 3-jet measurement function will 
make sure that all events where a singular limit 

— > is approached will be rejected. 

In the case of the one-loop contributions, the 
most complicated objects will be 5-point integrals 
with one off-shell external leg. Sector decompo- 
sition will lead to a result in terms of five inde- 
pendent scaled Mandelstam invariants y\. This 
result has to be calculated up to order e 2 , as it 
will be combined with the 1 — > 4 parton phase 
space where one parton can become unresolved, 
leading to 1/e 2 poles. This does in principle not 
constitute a problem, as the expansions to higher 
order in e, as well as 1 — > 4 parton phase space 
integrals, are well under control within sector de- 
composition. It is also possible to do parts of the 
loop integrations analytically to achieve a form 
which is suitable for subsequent sector decompo- 
sition ^Ej. However, these contributions have 
not yet been implemented into the Monte Carlo 
program. 

3.2. Real radiation at NNLO 

In [T^l, the method based on sector decompo- 
sition is for the first time applied to extract the 
poles appearing in massless 1 — > 5 particle inte- 
grals. The correctness of the results for the in- 
tegrals over the 1 — > 5 particle phase space can 



be checked by exploiting the fact that the sum 
over all cuts of a given (UV renormalised) dia- 
gram must be infrared finite. This is shown in 
Figurenfor a sample diagram: Summing over all 
cuts of this diagram and performing UV renor- 
malisation, we obtain the condition 

Ti->5 + z i T i^4 + z 2 Ti_» 3 + z 3 Ti_ 2 = finite , (3) 
where Ti_>i denotes the diagram with i cut lines. 
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Figure 1. Cancellation of IR divergences in the 
sum over all cuts of the renormalised graph 



The renormalisation constants Zi (in Feynman 
gauge) are given by |13ll8j 




The important new ingredient in eq. @ is the 
calculation of T\^. The sector decomposition 
method leads to [T%] 




where the numerical accuracy is 1%. The expres- 
sions entering eq. J3J for i < 5 combine to |TB] 



zi Ti^ 4 + z 2 Ti^ 3 + z 3 Ti^ 2 = 
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+ finite} . (6) 

We can see that the poles in © are exactly can- 
celled by the 5-parton contribution JSJ) within the 
numerical precision. 

3.3. Differential Monte Carlo program 

Although the sector decomposition approach is 
considered to be a "numerical method" , as the 
pole coefficients are only calculated numerically, 
the isolation of the poles is an algebraic proce- 
dure, leading to a set of finite functions for each 
pole coefficient as well as for the finite part. This 
feature allows the inclusion of any (infrared safe) 
measurement function, at the level of the final 
Monte Carlo program, which means that the sub- 
tractions and expansions in e do not have to be 
redone each time a different observable is con- 
sidered. However, in most cases of physical rele- 
vance, the measurement function is such that it 
would prevent certain poles from arising at all if 
it were included in the e-expansion. For exam- 
ple, poles associated with a 2-jet configuration, 
where 3 of the 5 final state particles become the- 
oretically unresolved, would be avoided ab initio 
by a 3-jet measurement function. As extra poles 
are rather costly in what concerns the creation 
of a large number of terms, it would therefore 
be desirable not to produce the terms associated 
with the subtraction of these poles at all in the 
e-expansion. In short, there is a certain dilemma 
between maximizing the flexibility to include any 
measurement function only at the stage of the fi- 
nal Monte Carlo program and limiting the size of 
the produced expressions. For e + e~ — > 3jets at 
NNLO, limiting this size is one of the most im- 
portant issues. Therefore, this dilemma has been 
solved by including some "preselection rules" in 
the e-expansion which reject configurations which 
will surely be 2-jet configurations. In the example 
shown here, this can be achieved by introducing a 
cut parameter y th - which must be smaller than 
any possible experimental resolution parameter 



corresponds to a 2-jet configuration. Note that 
this procedure does not introduce any dependence 
on y th in the matrix element, its only effect is 
to prevent the e-expansion subroutine from doing 
subtractions for S1345 — + 0. In this way, one can 
reduce the size of the expressions considerably 
without loosing any flexibility in what concerns 
the actual measurement function to be included 
later. In fact, the architecture of the program 
described in ^S] is the one of a partonic event 
generator, such that fully differential information 
about the final state is available. This architec- 
ture is independent of the actual matrix element, 
so the toy matrix element used in |18| and in this 
article can finally be replaced by a different mod- 
ule for the full matrix element, without destroying 
the property of providing the four-momenta of all 
final state particles. 

As an example, 3-, 4- and 5-jet rates using the 
JADE algorithm [221 are shown in Fig. [21 based 
on a toy matrix element built from the graphs 
shown in Fig. 
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Figure 2. 3-, 4- and 5-jet rates at order al for 
the toy matrix element 



4. Outlook 

The sector decomposition approach to the cal- 
culation of NNLO cross sections is a very powerful 
tool, especially in what concerns the double real 
radiation part, as it requires neither the manual 
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construction of subtraction terms, nor the factori- 
sation of the phase space and the analytic inte- 
gration of the subtraction terms in the singular 
limits. In |T2], the method is applied for the first 
time to a massless 1 — > 5 process. The code is 
constructed as a partonic event generator, which 
means that fully differential information about 
the final state particles is at one's disposal. 

A disadvantage of the sector decomposition ap- 
proach is given by the fact that it produces very 
large expressions, as in each decomposition step, 
the number of original functions increases. There- 
fore, CPU time is an issue for the treatment of 
processes with a large number of massless par- 
ticles in the final state, as for example e + e~ — > 
3jets at NNLO. However, the method sketched 
here relies on a division of the amplitude squared 
into different "topologies" corresponding to dif- 
ferent classes of denominator structures, such 
that the problem is naturally split into smaller 
subparts. Further, the size of the expressions can 
be reduced by including information about phys- 
ical limits already at the level of the e-expansion, 
without loosing any flexibility in what concerns 
the definition of observables at the Monte Carlo 
level. 

As the method is based on a universal al- 
gorithm acting on integration variables, it will 
surely see a number of interesting applications in 
the future. 
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